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Abstract: This paper addresses tlie distributed attitude synchronization problem of multiple spacecraft 
with unknown inertia matrices. Two distributed adaptive controllers are proposed for the cases with 
and without a virtual leader to which a time- varying reference attitude is assigned. The first controller 
achieves attitude synchronization for a group of spacecraft with a leaderless communication topology 
having a directed spanning tree. The second controller guarantees that all spacecraft track the reference 
attitude if the virtual leader has a directed path to all other spacecraft. Simulation examples are presented 
to illustrate the effectiveness of the results. 
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1 Introduction 

• In recent years, consensus and cooperative control of multi-vehicle systems have attracted com- 

tyj ■ pelhng attention from various scientific communities. A large body of theoretical advances has 

^ been reported, see [9l[8|ll H 113 1 114 1 [TO] s-nd references therein. In the aforementioned works, the 

^ ■ agent dynamics are restricted to be a single, double integrators or linear systems. The results 

K*" ' proposed in these papers become quite limited when dealing with the attitude synchronization 

■ problem of multiple spacecraft, which is more challenging than the consensus of vehicles with 

■ integrator dynamics, due to the nonlinearity of the attitude dynamics. 

^ ' Attitude control of a single rigid body has been extensively studied, e.g., in [22l EOl |2l |23] . A 

' leader- follower strategy is proposed in [21] for attitude synchronization of multiple spacecraft. 

Decentralized control laws using the behavioral approach are presented for attitude synchro- 
nization in [151 [7], where the communication topology among spacecraft is assumed to be a 
^ ' bidirectional ring. Adaptive consensus protocols are proposed in [3j for multiple manipulators 

with uncertain dynamics. Cooperative attitude control of multiple rigid bodies with a leader- 
follower communication topology is considered in [5j. In [1], contraction analysis theory is used 
to derive attitude synchronization strategies with global exponential convergence for a group 
of spacecraft with a bidirectional communication topology. Distributed control laws without 
velocity measurements are studied for attitude synchronization of multiple spacecraft in p!] by 
use of quaternion representation while in [17] by use of Modified Rodriguez Parameters (MRPs) 
for attitude representation. The results in [H [T7] are applicable to general undirected communi- 
cation topologies. The distributed attitude tracking problem is addressed in [16] for spacecraft 
whose information exchange graph can be simplified to a graph with only one node. 

Motivated by |171 [^[5]. this paper concerns the distributed adaptive attitude synchronization 
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problem of a group of spacecraft with unknown inertia matrices. The attitude dynamics are 
represented here by MRPs. Two distributed adaptive controllers are proposed for the cases 
with and without a virtual leader to which a time-varying reference attitude is assigned. The 
first controller achieves attitude synchronization for a group of spacecraft with a leaderless 
communication topology having a directed spanning tree. The second controller guarantees that 
all spacecraft track the reference attitude which is available to only a subset of the spacecraft, 
if the virtual leader has a directed path to all other spacecraft. Differing from the results 
given in [3] which applies only to a bidirectional ring communication topology, and those in [3] 
which requires the communication graph to be undirected, the communication topology among 
the spacecraft is relaxed to a general directed graph in this paper. The results obtained here 
generalize Theorem 4.1 in [17] to the case where the attitude dynamics are uncertain and the 
communication topology is either leaderless or leader- follower. It should be mentioned that all 
the results in this paper are applicable to robotic manipulators with dynamics represented by 
the Euler-Lagrange equation. 

The rest of this paper is organized as follows. The attitude dynamics and some useful results of 
the graph theory are introduced in Section 2. The distributed adaptive attitude synchronization 
problem for the cases without and with a reference attitude are studied in Sections 3 and 4, 
respectively. Simulation examples are presented to illustrate the theoretical results in Section 
5. Section 6 concludes the paper. 

The following notation will be used throughout the paper. R,"^" denotes the set of all re x n 
real matrices, /at represents the identity matrix of dimension N. 1 £ denotes the vector 
with all entries equal to one. Matrices, if not explicitly stated, are assumed to have compat- 
ible dimensions. ||M|| represents the induced 2-norm of matrix M € R"^™. For a vector 

-X3 X2 

X = [xi,X2,X3]-^ G R^, the cross-product operator is denoted by S{x) = X3 — xi 

—X2 xi 

diag(^i, • • • , An) represents a block-diagonal matrix with matrices Ai,i = 1, • • • , re, on its diag- 
onal. A® B denotes the Kronecker product of matrices A and B. 



2 Preliminaries and Problem Formulation 

This paper considers the attitude synchronization problem of a network of A'^ spacecraft. Mod- 
ified Rodriguez Parameters (MRPs) are used here to represent the attitude of the spacecraft 
with respect to the inertial frame. The MRP vector ai € R'^ for the i-th spacecraft is defined by 
cjj = eitan(^), where is the Euler axis and is the Euler angle [18]. The attitude dynamics 
of the 2-th spacecraft are given by [2U) 

JiLLli = -S{uJi)JiUJi + Ui, 

ai = G{ai)ui, i = 1,2,--- ,iV, 
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where G denotes the angular velocity in the body-fixed frame, Jj € R^^'^ is the inertia 
matrix, Uj € R^ is the control torque, and 

Following jl9l I23j . ([T]) can be rewritten as 

H*{ai)ai + C*{ai,a^)ai = G-^iai)ui, i = 1, 2, • • • , TV, (2) 

where 

H:{ai)^G'^{ai)J,G-Hai), 
C:{a„&i) ^ G-^{a,)J^G-\ai)G{a,)G-\a,) - G-^ {a,)S{JiG-\ai)&i)G-\ai). 

It is assumed that the inertia matrices Jj, i = 1,2, •• • ,N, are unknown constant positive- 
definite matrices. Under this assumption, the attitude dynamics ([2]) has the following properties: 
Property 1. Matrix H*{ai) is symmetric and positive definite. 
Property 2. Matrix H*{ai) — 2C*{ai,&i) is skew symmetric, i.e., 

x^{H*{ai) -2G*{ai,ai))x = 0, \/xeK^. 

Property 3. The attitude dynamics ([2]) satisfies the following linear parameterization condi- 
tion: 

H*{(Ti)y + C*{ai,ai)y = Y{ai,ai,y,y) 9^, Vy G R^, 
where Y G R'^^^ is called the regression matrix and 

1 T 



'11 "^12 "^13 "^22 "^23 "^33 



(3) 



is the unknown parameter vector with Jjf^ being the (j, k)-th entry of the inertia matrix Jj in 

The communication topology among spacecraft is represented by a directed graph Q consisting 
of a node set V = {1, 2, • • • , A^} and an edge set £ C V x V. The node i denotes the z-th spacecraft. 
An edge means that spacecraft j can obtain the attitude information of spacecraft i, but 
not conversely. For an edge in the directed graph, i is the parent node, j is the child node, 
and j is neighboring to i. A graph with the property that G £ implies (j, i) G f is said to 

be undirected. A path on Q from node ii to node ii is a sequence of ordered edges of the form 
{ik, ik+i), k = 1, - ■ ■ ,1 — 1. A directed graph contains a directed spanning tree if there exists a 
node called the root, which has no parent, such that there exists a directed path from this node 
to every other node. 

The adjacency matrix A G R^^^ of graph Q is defined as an = 0, and aij > if G £ 
but otherwise. The Laplacian matrix C G R^^^ is defined as Cu = Ylj^ti^iji ^ij — 
i ^ j- Given a matrix R = {rij)pxp, the graph of R is the directed graph with p nodes such that 
there is an edge in the graph from node j to node i if and only if rij 7^ [6]. 

Lemma 1 [13]. Zero is an eigenvalue of C with 1 as the corresponding right eigenvector and 
all the nonzero eigenvalues have positive real parts. Furthermore, zero is a simple eigenvalue of 
C if and only if the graph has a directed spanning tree. 



3 



3 Distributed Adaptive Attitude Synchronization 



This section considers the distributed adaptive attitude synchronization problem of ([T]) whose 
communication topology is represented by a leaderless directed graph Q. A graph is leader- 
less, if each node in this graph has at least one parent. Before moving forward, the attitude 
synchronization problem is first defined. 

Definition 1. The distributed adaptive attitude synchronization problem is said to be solved, 
if the control laws Ui, i = 1,2, ■■ ■ ,N, are designed by using only local information of neighboring 
spacecraft such that the attitudes of ([T|) satisfy limf_j.oo — crj\\ = 0, linit^oo — o'jW = 0, 
i,j = 1,2,- •• ,N. 

At each time instant, the attitude information of other spacecraft available to spacecraft i is 
given by 



CTi — 



2^ 1=1 



(4) 



where A = {aij)is[xN is the adjacency matrix of the communication graph Q. 

To quantify whether the attitude synchronization is achieved or not, a synchronization error 
ei{t) is defined as follows: 

,iV. (5) 



ei = ai-af, i = l,2. 



In addition, further define a filtered synchronization error Si{t) as 

Si = ei + Aid, i = 1, 2, • • • ,N, 

where Aj G R^^^ 



(6) 



i = 1,2, ■ ■ ■ ,N, are constant positive-definite matrices. 
In light of (l2|), ([5]), and dGj), it can be obtained that vector Si evolves according to the 
following dynamics: 



H*{ai)s, + C*{ai,a,)s, = C {a,)ui- H*{ai){af - A,e,) -C*{ai,a,){&f - AiCi), i = l,2,--- ,N. 

(7) 

By Property 3, the right side of the above equation can be written into a linear combination of 
the inertia vector 9i, which is defined in ([3]). To the end, introduce a linear operator L{a) for 
vector a = [01,02,03]"^ as 

oi 02 as 
L{a) = oi 02 03 
oi a2 (13 

and a linear operator F{x,v) for vectors x = [xi,X2,X2]'^ , v = [^1,^2,^3]'^ as 



F{x, v) 



X1V3 -X1V2 X2V3 -X2V2 + X3V3, -X3W2 

-X1V3 -X2V3 xivi - X3V3 X2V1 X3V1 

X1V2 -xivi + X2V2 X3V2 -X2V1 -X3V1 



It can be verified that operators L{s) and F{x,v) satisfy 

Jia = L{a)9i, S{Jix)v = F{x,v)9i. 



(8) 
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Therefore, by using ([5]), ([7]) can be written as 

H*{ai)si + C*{ai,&i)si = G-^{ai)ui -Yi{ai,ai,&f - Aiei,af - AiCi) Oi, i = 1, 2, • • • , iV, (9) 
where 

Yi{a^,&^,af - Aieuaf - Kei) = G-^((J.) [L{G-\cji){af - A,e,)) + L{G-\ai){af - A,e,)) 

+ F{G~\ai)ai,G-\a^){df - Kei))) 

Since the inertia parameter 9i is unknown, its estimate 6i is used instead to construct the 
controller Ui to spacecraft i as follows: 

Ui = G'^{ai){Yi{ai,ai,af - Kiei,df - kiei)ei- KiS^ , i = 1, 2, • • • , iV, (10) 

where Ki G R^^^ i = 1, 2, • • • , A^, are positive definite. The parameter estimate vector 9i is 
generated by the following adaptive updating law: 

Oi = -TiY^{ai, &i,af - Kei^af - Kiei)si, i = 1, 2, • • • , TV, (11) 

where Fj G R^^^, i = 1, 2, • • • , A^, are positive-definite diagonal matrices. 

Define the estimation errors 9i = 6i — 9i, i = 1,2, ■ ■ ■ , N . Then, substituting (jlOp . (jlip into 
dS]) gives 

H*{ai)si + G*{ai,&i)si = -Yi{ai,ai,&f - Aiei,af - Aiei)6i- KiSi, i = l,2,--- ,N. (12) 

Theorem 1. If the leaderless communication graph Q has a directed spanning tree, then 
distributed controllers (jlOp and adaptive updating laws (jlip solves the attitude synchronization 
problem for ([1]). 

Proof. Let e{t) = [ej , ■ ■ ■ ,e^]"^ and a{t) = [af , ■ ■ ■ ,ajf]'^. Then, ^ can be written as 

e = MCa, (13) 

where £ is the Laplacian matrix associated with graph Q, and M = diag(^^;^ ay, • • • , Z^jLi ciNj) 
Since graph Q is leaderless and has a directed spanning tree, one obtains 1) there exists at least 
one nonzero entry for each row of the adjacency matrix A, thus matrix M is positive definite; 
2) the Laplacian matrix £ has a simple eigenvalue with 1 as the corresponding eigenvector, 
and the other eigenvalues have positive real parts. Then, it follows from (|13p that e = if and 
only if ui = (72 = • • • = ctat. That is, the attitude synchronization problem is solved if and only 
if e{t) 0, e{t) 0, as t — ;> oo. 

Consider the following Lyapunov function candidate 

/TV N \ 

V{t) = 2 E 'lH:{a,)s, + ^ 9[Tr% . (14) 

\i=l i=l / 
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By (|12p . pip , and Property 2, the time derivative of V{t) is 

N N 
TV 

= ^s^ [H*{ai)si - C*{ai, &i) - Yi{ai, &i, &f - Kei^af - Aiei)0i - KiSi^ 

i=l 

N 

(15) 

i=l 
N 

i=l 

N 

i=l 

Let S = {{(Ji,Si,ei)\V = 0}. Note that F = imphes that Si = 0, i = 1, 2, • • • , iV. By 
LaSalle's invariance principle [19], it follows from that Sj — )• 0, i = 1, 2, • • • , A^, as t — )• cx), which 
by dHl) in turn shows that e{t) — )• 0, e{t) — )• 0, as t — )• oo, i.e, attitude synchronization is achieved. 
■ 

Remark 1. It should be noted that the controller (jlOp to spacecraft i depends only on 
its own attitude vectors Ui, &{, a-i, and the attitudes of its neighboring spacecraft, therefore 
is distributed. Theorem 1 gives a sufficient condition for achieving attitude synchronization. 
However, it is generally quite hard to expressly derive the final synchronized attitude value, 
which depends on the initial attitudes of the N spacecraft, matrices Aj, iCj, i = 1, 2, • • • , A^, and 
the communication topology Q. 

4 Distributed Adaptive Attitude Tracking 

Different from the leaderless communication graph discussed as in the above section, the space- 
craft' attitudes may be desired to follow a given time- varying reference attitude u'' in certain 
circumstance. It is supposed that a"^ is available to only a subgroup of the spacecraft, other- 
wise cooperation between neighboring spacecraft by exchanging attitude information become 
not so necessary. Assume that o"'", cr^', and are all bounded. For this case, the attitude 
synchronization problem is called attitude tracking problem in [H [T7] , which is formulated as 
follows. 

Definition 2. The distributed adaptive attitude tracking problem is said to be solved, if the 
local control laws Uj, i = 1,2,- •• ,A^, are designed for ([T]) such that limj__j.oo — a'^\\ = 0, 
limt^oo - fj^'ll = 0, i = 1,2, • • • ,7V. 

Take the reference attitude a''' as the attitude of a virtual leader, labeled as spacecraft + 1. 
Since the virtual leader does not obtain any information from the spacecraft, the communi- 
cation topology among these A^ + 1 spacecraft (the A^ spacecraft and the virtual leader) is in 
the leader-follower form. 
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Assume that the communication topology among the N spacecraft is still denoted by Q. The 
attitude information of other spacecraft available to spacecraft i is given by 

4 i = 1.2,...,iV, (16) 

where A = {aij)iyxN is the adjacent matrix of G, aj(7V4-i) > 0, i = 1, • • • , A^, if spacecraft i 
has access to the virtual leader and aj(7v+i) = otherwise. Similar to the above section, a 
synchronization error efi{t) is defined as follows: 

efi = afi-a%, i = l,2,---,N. (17) 

Correspondingly, the filtered synchronization errors s/j is defined as 

Sfi = efi + AiCfi, i = l,2,---,N, (18) 

where matrices Aj E R^^^, i = 1,2, - ■ ■ ,N, are positive definite. By ([2]), p^ . (fT7|) . and Property 
3, vector Sfi satisfies the following dynamics: 

Ht{(yi)sfi + C*{ai,ai)sfi = G~'^{ai)ui - Yi{ai,&i,a'}i - Aiefi,a'ji - Aje/j) 9i, i = 1,2, - ■ ■ ,N. 

(19) 

The distributed controllers Ui, i = 1,2, ■■ ■ , N, to the N spacecraft are proposed as 

Ui = {Yi{ai,ai,dji-Kiefi,aji-Kiefi)6f^-KiSf^ , i = l,2,--- ,N, (20) 

where Ki € R^^^, i = 1,2, - ■ ■ ,N, are positive definite and 9fi is the estimate of 9i, generated 
by the following adaptive updating law: 

9fi = -fiY^ {audi, aji- lief i,aji- lief i)sfi, i = l,2,---,N, (21) 

where Fj € R^^^, z = 1, 2, • • • , A^, are positive-definite diagonal matrices. 

Theorem 2. Denote by ^at+i the directed graph of matrix A]\jj^i = g q ' ^^^^^'^ ^ ~ 

[ai(jv+i)) ■ ■ ■ 1 o.N{N+i)Y' ■ If graph Gn+i has a spanning tree with node A + 1 as the root, then 
distributed controllers (j2U|) and adaptive updating laws (|21|) solve the attitude tracking problem 
for the A spacecraft in ([1]). 

Proof. Let gn+i = o"'", a{t) = [cjf , • • • , cj^, a^+J"^, and ef{t) = [eji,--- ,ejjv,0]'^. Then, 
(jl7p can be written as 

ef = MCa, (22) 

where M = diag(^^j^ aij, • • • , X^^i a^j, 1), and £ = {Cij){N+i)x{N+i) is the Laplacian matrix 
of graph §, defined as Cu = Ylj'=i ^ij^ ^ij = -^ij U / 0) G {1, 2, • • • , A}, and C(^N+i)j = 0, 
Vj € {1,2, • • • , A^ + 1}. Under the assumption of the theorem, matrix M is positive definite 
and is a simple eigenvalue of matrix C with 1 as the corresponding eigenvector, implying that 
Cf = if and only if ai = (T2 = • • • = (Tn = cr^ ■ That is, the distributed attitude tracking 
problem is solved if and only if e/(t) — )• 0, e/(t) — )• 0, as t — )• oo. 
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Consider the following Lyapunov function candidate 

^ 1 / ^ ^ ^ \ 

\i=l i=l / 

The time derivative of V{t) can be obtained as 

N N 
1=1 

N 



(23) 



i=l 



(24) 



i=l 



Since V{t) > 0, y(t) < 0, V{t) remains bounded, which by (i23]l implies that Sfi, 9fi, i = 
1, 2, • • • ,N, are bounded. By using standard signal chasing arguments, it is easy to show that 
s fi, i = 1,2, N , are bounded. Therefore, V{t) = is bounded. In light of 

Barbalat's lemma [19], V{t) ^ 0, as t — )• oo. Thus, Si — )■ 0, i = 1, 2, • • • , N , as t ^ oo, which by 
(jl7p shows that e/(i) — )■ 0, e/(t) — > 0, as t — >• oo, i.e, the distributed attitude tracking problem 
is solved. ■ 
Remark 2. Theorems 1 and 2 present sufficient conditions for the adaptive attitude syn- 
chronization of multiple spacecraft having general directed communications in the presence of 
unknown inertia matrices, for both the cases with and without a reference attitude. By contrast, 
the results given in [2] are applicable only to a bidirectional ring communication topology, and 
the result in [3j requires the communication graph to be undirected. Theorems 1 and 2 gen- 
eralize Theorem 4.1 in |17) to the case where the attitude dynamics ([T]) are uncertain and the 
communication topology is either leaderless or leader-follower. It should be noted that different 
from Theorem 1, Barbalat's lemma is utilized to derive Theorem 2, due to the fact that the 
time-varying reference attitude a^' may render the closed-loop system nonautonomous. 



5 Numerical Examples 

In this section, the effectiveness of the proposed control laws is illustrated through numerical 
simulations. 




Figure 1: The communication topology. 
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Consider a group of six spacecraft, whose inertia matrices are shown in Table 1 [T5J. The 
initial attitude states cjj(O), Wi(0), i = 1,2, •■■ ,6, are chosen randomly. The communication 
topology is given by Fig. 1, so the corresponding adjacency matrix is 



A = 












1 


1 


1 


1 

















1 


1 














1 


























1 




















1 






For ([TO]) and 1^, take matrices = = I3, Ki = Ki = 3/3, Pj = f j = Sle, i = 1, 2, • • • ,6. The 
parameter estimates 9i and 9fi are initialized to be zero, i.e., 0i(O) = 0, Ofi{0) = 0, « = 1, 2, • • • ,6, 
in (fTTj) . ([2T]) . For simplicity, let the reference attitude a'"" = [0.1,0.3,0.5]-'". Suppose that is 
available only to spacecraft 1. In this case, an = 1, 0^7 = 0, i = 2, ■ ■ ■ ,6. 



Jl 


[1,0.1, 0.1; 0.1, 0.1, 0.1; 0.1, 0.1, 0.9] kgm^ 


J2 


[1.5, 0.2, 0.3; 0.2, 0.9, 0.4; 0.3, 0.4, 2.0] kgm^ 


h 


[0.8, 0.1, 0.2; 0.1, 0.7, 0.3; 0.2, 0.3, 1.1] kgm^ 




[1.2, 0.3, 0.7; 0.3, 0.9, 0.2; 0.7, 0.2, 1.4] kgm^ 




[0.9, 0.15, 0.3; 0.15, 1.2, 0.4; 0.3, 0.4, 1.2] kgm^ 




[1.1, 0.35, 0.45; 0.35, 1.0, 0.5; 0.45, 0.5, 1.3] kgm^ 


Table 1: Spacecraft specifications. 



Figs. 2(a), 2(b), and 2(c) depict, respectively, the attitudes, angular velocities, and control 
torques of the six spacecraft with and pT|) . from which it can be observed that attitude 
synchronization is indeed achieved. The parameter estimates Oi, i = 1,2,-- - ,6, are shown in 
Fig. 3. Figs. 4(a), 4(b), 4(c), and 5 depict, respectively, the attitudes, angular velocities, control 
torques, and the parameter estimates Ofi, i = 1,2, ••• ,6, of the six spacecraft with (f20jl and 

(EH). 



6 Conclusions 

This paper has addressed the distributed adaptive attitude synchronization problem of a group 
of spacecraft with unknown inertia matrices. Two distributed adaptive controllers have been 
proposed for the cases with and without a virtual leader to which a time- varying reference atti- 
tude is assigned. The first controller achieves attitude synchronization for a group of spacecraft 
with a leaderless communication topology having a directed spanning tree. The second controller 
guarantees that all spacecraft track the reference attitude if the virtual leader has a directed 
path to all other spacecraft. This paper has extended some existing results in the literature. 
An interesting topic for future research is the distributed adaptive attitude synchronization of 
multiple spacecraft without velocity measurements. 
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Figure 2: Attitudes, angular velocities, and control torques of ([T|) with controller (|1U|) and 
adaptive law (llip . 
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(a)0„i = l,2,3 (b) e',, i = 4,5,6 

Figure 3: Parameter estimates of (jlip . 
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(c) Control torques 



Figure 4: Attitudes, angular velocity, and control torques of ([T]) with controller (|2Up and adaptive 
law (EB. 
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Figure 5: Parameter estimates of (j2ip . 
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